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1.0 INTRODUCTION 


1.1 OBJECTIVE 

The original objective of this program was to perform studies that would develop 
enhanced contingency power capability for turboshaft engines in tiltrotor applications. 
The intent was that the result of these studies would yield a design definition or 
methodology that could lead to a reduction in both engine size and the cost required to 
hover safely on one engine thereby improving aircraft economics. The measure to be 
used in these studies was: Reduced aircraft operating cost versus conventional 
turboshaft. 

1.2 PHASE I STUDY RESULTS 

Phase i of this program was completed and reported on NASA Contract Number NAS3- 
26617, Task Order #48. In phase I, a study plan was put together as follows: 

1 . ) Establish baseline engine and aircraft 

2. ) Identify and screen contingency power and enhancing concepts 

3. ) Select best concepts, define, and evaluate concepts versus baseline 

4. ) Identify key technologies and development plan. 

From this plan, the following two concepts were selected for further definition and 
evaluation: 

1 . ) Throttle push with cooled power turbine 

2. ) Modulated high pressure and power turbine cooling air. 

The results of these studies were as follows: 

« Modulated cooling and throttle push concepts yield biggest payoff and have highest 
probability of success. 

- 3 to 4% improvement in DOC => +30% to +80% in operating profits. 

• For the studied range of compressor pressure ratios (1 .25 to 1 .4), both concepts 
have similar payoff. 

- Best concept for application depends on specific engine configuration and 
mission requirements. 

The study was performed using standard lifing methodologies which can’t accurately 
model short duration material properties. From this, it was concluded that a key 
technology for most contingency power concepts is the development of a life prediction 
rnethodoiogy for high power, short duration ratings. It is believed that better modeling 
could yield further DOC improvements and iower development program cost and risk 
and would allow existing engines to be rerated for higher contingency power. This 
formed the basis for the phase II of this program which is reported on here. 
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1.3 PHASE II PROGRAM 


Phase I! of this program was aimed at improving materials knowledge and fifing 
methodologies. The goals of the program were as follows: 

1 . ) Conduct testing to characterize materia! behavior at these conditions 

2. ) Develop a computational life prediction methodology integrating the obtained 

material data. 

1 .3.1 Phase II - Material Test Program 
Material Description 
Utilizes two advanced alloys: 

1 . ) Rene' 142 Directionally Solidified Alloy (DS) - prime HPTB material used in many 

GEAE applications including CF6, CFM56, CF34, CT7, and GE90 

2. ) Rene’ N6 Single Crystal Alloy (SX) - 3rd generation SX alloy also in production. 

Both alloys are backed by a large mechanical and physical property data base that is 
practically truncated at the contingency power setting condition. 

Testing Philosophy 

1 . ) Follow test approach of Arnold, Castelli and Saleeb to generate GViPS materials 

constants, which was a challenge for the commercial test vendors 

2. ) Test in the longitudinal orientation - principal stress direction 

3. ) Use time and temperature frame consistent with OEI event 

4. ) Use closed loop computer controlled test machines “record all the data” 

5. ) Simulate an OEI event in test environment 

Figure 1-1 shows the Phase II test matrix and Figure 1-2 shows the validation tests. 
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Stress Relaxation/Tensile 
Temperature 160Q°F-23Q0°F 
Strain Rate 1.0E-06 - I.OOE-OI/s 
1 22 tests 



Constant Stress Creep 

Temperature 1 600°F-2300°F 
Times 1-24 hrs plus some longer 
50 tests 



Stepped Stress Creep 

Temperature 1600°F-2300°F 
22 total 




Constant Stress Creep 
Temperature 1600°F-2300°F 
Times 1- 24 hrs plus some longer 
50 tests 



Stepped Stress Creep 
Temperature 1600°F-2300°F 
22 total 



Figure 1-1. Phase !! Test Matrix. 
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SPLCF 


Temperature 1800°F-2200°F 
1 6 tests 



co 


Time 


Cyclic Temperature Creep - Simulates the 0E1 Event 

Temperature 2000°F - 2200°F 
1 2 tests 


Temp. 


2000F 


2250F/2 min. 


2000F 


Temp. 




2250F/30 s. 



2150F/2 min. 

-100-200 hr. failure 


Time 


Time 


Figure 1-2. Validation Tests. 

1 .3.2 Life Prediction Methodology 

Locations of typical turbine blade life limiting mechanisms are shown in Figure 1-3. 
Currently, turbine blades are iifed for Stress Rupture (creep) or for Low Cycle Fatigue 
(LCF). Tensile yield is evaluated for engine over-temperature certification. Turbine 
blades are assumed “damaged” after a single over-temperature event. No attempt is 
made to quantify over-temperature damage with respect to stress rupture or LCF life 
limits, which is a practice that introduces an unnecessary amount of conservatism into 
the blade design. In addition, the current blade lifing approach de-couples inelastic 
strain as outlined below: 


Current Approach 

^total 

8 


elastic elastic* creep 

= 8 +8 + 8 _» 


Time independent analysis assumes: 


'Time Independent 
-Time Dependent 


8 


total 


Inelastic Strain = Plastic Strain (Tensile Yield) 

Design requirement: survive an 
over-temperature event without 


elastic plastic 

8 "he 


tensile yield leading to component 
failure (short time). 
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Time dependent analysis assumes: 


Inelastic Strain = Creep Strain (Stress Rupture) 

Design requirement: endure 
application mission life without 

C3TS6P 

S tertiary creep (onset of voiding) 

in component (long time). 


„ total 


= s 


elastic 



Tip oxidation 
(time at temp, 
and cycles) 

Low Cycle 
Fatigue (cycles) 


Tensile yield 
(over- temp.) 


Stress Rupture 
(time at temp.) 


Figure 1.3. Turbine Blade Life Limiting Mechanisms. 


The new life prediction methodology developed in phase II integrates the material data 
from the test program with a unified approach that couples inelastic strains as described 
below: 


Unified Approach 


£ 


total 


elastic plastic creep 

8 +8 +8 


8 


total 


£ 


elastic 


+ 



£ 


inelastic 


Modeling material response to over-temperature events will allow prediction of impact on 
residual stress rupture life. 

» Several unified constitutive models are available. Models evaluated in this study are: 

- Generalized Viscoplasticity with Potential Structure (GVIPS), Arnold and Saleeb 

- Single Crystal, Stouffer, Dame, and Sheh 
Both Models Developed Under NASA Funding 

Details of the test program and the life prediction methodology are discussed in the 
following sections. 
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2.0 MECHANICAL TEST PROGRAM 


2.1 INTRODUCTION 

The ability to design advanced turbine components critically hinges on the description of real 
material behavior. Since turbine blades and nozzles experience wide spectrums of temperature 
and loading, predicting component response requires design methodologies which vary from the 
classical approaches of independent plasticity and creep. Additionally, in order to attain 
optimum strength these material descriptions must address the anisotropy inherent in modern 
Ni-base superalloy components, since these usually have either a directionally solidified (DS) or 
single crystal (SX) microstructure. 

The problems cited have resulted in considerable interest in the “state variable” approach to 
material constitutive behavior. This approach has developed from a synthesis of modeling, 
physical metallurgy and experiment over the last two decades 2 . The historical development of 
these state variable models can be traced in many references 3 ' 14 . A multitude of models exist 
and several state variable-based models have been applied to Ni-base superalloys 15 ' 26 ; 
however, two models are of particular interest to the OEI Tiltrotor initiative. These are the 
Generalized Viscoplasticity and Potential Structure (GV1PS) approach of Arnold, Saleeb and 
Casteili 22 ' 26 and the Single Crystal models of Stouffer et al 17 ' 19 . Potentially, the GVIPS mode! 
can be applied to both DS and SX materials while the Single Crystal model is unique to the 
single crystal deformation restraints. Particular attention has been made in the GVIPS 
approach to describe unique material mechanical tests, which can directly relate to the materia! 
parameters needed in the model 22 ' 26 . 

2.2 APPROACH 

The mechanical test program was based largely on the GVIPS testing rational, using both an 
advanced directionally solidified alloy and its single crystal counterpart. The alloys tested were 
directionally solidified Rene 142 and single crystal Rene N6. Rene 142 is described in US 
patent 5,173,255 27 while the RN6 description is embodied in US patents 5,270,1 23 28 and 
5,455, 120 29 . 

All tests were performed from the same production qualified master heat for each alloy. 

Castings were provided by PCC’s Development Facility in Minerva, Ohio. Both 0.5 inch 
diameter rods and 3/8 inch thick plate were cast. Single crystal primary orientations were held 
to within 15 degrees. 

The specific test goals were as follows: 

s Provide a basis for the development of the model parameters for the GVIPS and Single 
Crystal models and test the validity of these treatments. 

« Demonstrate the utility of the Casteili, Arnold and Saleeb 24 test methodologies in Test . 
Laboratories with advanced testing capabilities while extending these tests to the high 
temperature regime of OEi operation 

• Design a test to simulate an OEI event. 
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Two GE approved Testing Laboratories were utilized in the program. These were Material 
Characterization Laboratory (MCL) in Scotia, NY, and Metcut Research Incorporated in 
Cincinnati, OH; both are routinely utilized by General Electric for critical mechanical tests. 

2.3 TENSILE/RELAXATION TEST PROCEDURE 

Three types of tensile/stress relaxation tests were performed at MCL including: 

1 A.) Tensile test strained to 2% strain followed by 2 hours relaxation. The 1 A test matrix 

included testing at 1600°F, 1800°F, 2000°F, 2100°F, 2200°F, 2250°F and 2300°F at strain 
rates of 1 .OE-06/s, 1 .OE-05/s, 1 .OE-03/s and 1 .OE-OI/s. Rene 1 42’s upper test limit was 
2250°F. 

I B. ) Tensile test strained to 1 % strain followed by 2 hours relaxation, then strained to 2% 

followed by 2 hours relaxation. The 1 B test matrix included testing at 1 800°F, 2000°F, and 
2200°F at strain rates of 1 ,0E-05/s, 1 .OE-03/s and 1 .OE-OI/s. Rene 1 42’s upper test limit 
was 2250°F. 

I C. ) Tensile test strained to 0.5% strain followed by 2 hours relaxation, then strained to 1 % 

followed by 2 hours relaxation, then strained to 2% followed by 2 hours relaxation. The 1 C 
test matrix included testing at 1600°F, 1800°F, 2000°F, 2100°F and 2200°F at strain rates 
of 1 .0E-05/S, 1 .0E-03/S and 1 .OE-OI/s. 

An example of the three tensile relaxation test types is given in Figure 2-1 for R142 at 1 800°F 
and 1 .OE-05/s strain rate. 


Tensile/ Relaxation Test Types, R142 Alloy, I.OE-05/s 



Strain (in/in) 

Figure 2-1. Rene 142 Tensile Relaxation Behavior At 1800°F And I.OE-05/s Strain Rate. 
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MCL’s test specimen geometry was a standard threaded 0.200 inch diameter specimen with a 
0.750 inch smooth gage section. These were machined from the cast rod stock. The 
specimens were bare, since both R142 and RN6 have excellent oxidation resistance. MCL 
used standard MTS servo-hydraulic LCF type test stands at strain rates of 1 .0E-03 and greater. 
Low-pressure surface contact type extensometry was used to measure strain and provide the 
control signal for constant strain rate control. The slow strain rate tests were performed on MCL 
custom designed electro-mechanical test machines. These use a stepper motor to turn a ball 
screw which drives an actuator up and down and can be programmed for stroke, load or strain 
control. The same extensometer set-up is used with these machines. 

Metcut performed a similar test matrix except, generally without relaxation, and about 50% of 
the tests were run to failure. Metcut utilized a standard threaded 0.1 60 inch diameter specimen 
with a 0.750 inch gage section. These specimens were machined from the cast plate stock and 
were coated with a Pt-AI diffusion coating. Note that the primary dendrite spacings and 
microstructures for plate and rod were comparable. Metcut’s tests were performed on a 
standard Instron screw type machine using Instron Series 9 software. The extensometer set-up 
was similar to MCL’s and both vendors maintained a maximum temperature delta along the 
gage of 2°F through 2100°F and 3°F above this temperature. 

2.4 TENSILE/RELAXATION TEST RESULTS 

In general there was excellent agreement between both vendors as Figure 2-2 confirms. This is 
a plot of the 0.2% flow stress for Rene N6 over the entire temperature strain rate space. 



Figure 2-2. Rene N6 0.2% Yield Strength Behavior. Data From Both Test Vendors Is Plotted. 

Each vendor encountered problems. At MCL, data recording and control issues were 
encountered, especially at the 1 .OE-01/s strain rate. Here it was discovered that the control 
computer required a CPU with equivalent speed to at least a 166 MHZ Pentium chip. MCL 
logged data at a 100 points per second rate while Metcut utilized 50 points per second. 
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Both vendors ran into load train resonance at critical temperature - strain rate thresholds. 
Metcut, with the smaller test specimen, encountered vibration, set in motion by the machine’s 
acceleration, at 1 800°F and 1 .0E-1/s. The Metcut fix was to switch to “position” control, which 
afforded a slower ramp up rate at the start of the test. In position control, an “end cutter” or 
location sensor on the crosshead controls crosshead speed. Strain rate is calculated using the 
effective gage length. The strain is still measured by the extensometer. Figure 2-3 shows a 
comparison between extensometer strain rate control and position strain rate control in the 
resonance regime. 


R142 2000F Tensile Tests, Extensometer vs. Position Control 



— a — Modulus 
— 0 — # 15-1 Position 

#15-8, Extensometer 
o #15-1, Strain Rate (1/s) 
— 9 — #15-8, Strain Rate (1/s) 


0 ) 

« 

£E 


c 

GO 


Figure 2-3. Two Tensile Tests From Metcut Comparing Extensometer Versus Position Strain 
Rate Control In Rene 142 At 2000° F And Nominal 1.0 E-01/S Strain Rate. 

For MCL this problem was encountered at 2000°F and 1 .OE-03/s. At MCL a switch to the 
electro-mechanical machine helped to reduce but not eliminate the problem. Unfortunately, the 
electro-mechanical apparatus has an upper strain rate limit of 1.0 E-03/s, so the problem 
persisted at the highest strain rate conditions. Here, data smoothing routines helped in the 
more moderate resonance cases. It is believed that an appropriate strain rate acceleration 
schedule or some system damping design would have been of benefit. 

Another nuisance encountered at high temperatures was a cyclic “noise” signal picked up while 
measuring the low loads encountered at elevated temperatures. As such it was most 
pronounced during stress relaxation. Figure 2-4 shows a typical event, where the top curve is 
the data and the bottom curve is just the isolated noise. These events were intermittent and 
random in nature. No noise source was ever identified, and often a simple sine curve provided 
an excellent fit to the phenomenon. 
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R142 Test #34 Relaxation "Noise", 2200F, IE-05 per s. 



o Load * 

cal. £ 

y2-y1 

A y2-yi 


Relaxation Time (hrs) 

Figure 2-4. Cyclic Noise Occasionally Noted During Tensile Relaxation Test At Low Loads. 


2.5 CREEP TESTS 


MCL’s electro-mechanical test machines offered constant stress creep capability with controlled 
strain rate loading schedules. Data was accumulated from the true test start, and all plastic 
deformation was recorded. This test apparatus was well suited for stress change testing. An 
example of constant stress creep response of R142 at 1800°F is given in Figure 2-5. Also 
included in this figure is a cyclic stress change test. 

Metcut performed standard constant load creep tests per ASTM El 39 using the Vertex data 
control and logging system. Loading curves are automatically calculated and plastic strain 
recorded. Strain measurements were performed with transducers with a tolerance less than 
0 0001 inch. Since the true stress increases with extension in the constant load tests, these 
prove to have shorter lives as Figure 2-6 shows. Overall, the MCL test technique proved more 
reliable and provided less creep scatter. 


The following test matrix represents the final valid data content suitable for model analysis. 
Defective tests (with incompletely recorded data, extensometer/transducer slippage and 
equipment failure) have been set aside. 


Temperature 

Constant Load R1 42 
Constant Stress R1 42 
Constant Load RN6 
Constant Stress RN6 


1600°F 

12 

8 

9 

7 


1 800°F 

11 

9 

8 

6 


200Q°F 2100°F 2200°F 

Number Of Test Points 
10 12 9 

7 3 8 

9 8 9 

7 3 4 


2250 °F 2300°F 

9 8 

3 3 

9 7 

3 3 
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R142 Constant Stress Creep and Cyclic Stress Creep - 1800F 



Figure 2-5. Rene 142 Constant Stress Creep Data At 1800°F. 


R142 2200F Rupture Summary, Constant Load and Stress Tests 



0 — ' 

0.1 1 10 100 

Time (hrs) 


Figure 2-6. Rene 142 Creep Rupture Behavior At 2200°F. 

2.6 MIXED CREEP-PLASTICITY TESTS 

Since the state variable approach treats classical plasticity and creep as inelastic strain, a 
“mixed mode” test type can prove useful in validation of the procedure. This means that starting 
a test either first in a tensile mode then switching to creep or vice versa becomes a bench mark 
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test. One test of each type was completed at 1600°F, 1 800°F, 2000°F and 2200°F. An 
example of this type of test is given in Figure 2-7. 



Time (hrs) 

Figure 2-7. Rene N6 “Mixed Mode” Test Starting In Tensile Deformation, Switching To 
Creep, Then Ending In Stress Relaxation. 

2.7 SIMULATED OEI EVENT 

A special test was designed to simulate the brief OEI over temperature event. This test was set 
up on an MCL electro-mechanical test stand fitted with a custom designed igniter furnace. Care 
was taken to control temperature during the brief over temperature “spikes”. Basically the test 
was started in the normal constant stress mode. Just prior to the over temperature, the mode is 
changed to load control. After the temperature spike is complete and once thermal stabilization 
has occurred, the mode is switched back to constant stress. Multiple spikes at different 
temperature levels can be programmed and executed. Figure 2-8 shows a complete test 
performed on R142 where four 200°F over temperature simulation “spikes” were performed and 
after the last spike, the temperature was raised 1 00°F. Figure 2-9 shows a double spike event 
on an expanded scale. In this figure, the thermal expansion recovery is apparent. A total of 
eleven tests were performed with 9 on Rene 142 and 2 on Rene N6. 

2.8 RESULTS AND CONCLUSIONS 

The test types described by Castelli, Arnold and Saleeb 24 can be performed at commercial test 
sources which have state-of-the-art test technology. Several key points were learned during 
this effort including; 

* Control software expertise is paramount 
9 Resonance problems at high strain rate need to be addressed 
9 Noise and signal to noise limitations at high temperatures are always present 
9 Control computers must be fast and bug free. 
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Renel 42 Cyclic Simulated Over Temperature Event 



Figure 2-8. Rene 142 Simulated OEl Test 


Renel 42 Cyclic Simulated Over Temperature Event 
"Double Spike 1 ' 



Strain 

□ Temparatura 


Time (Min.) 


Figure 2-9. Consecutive 2 And 3 Minute Temperature Spikes. The Thermal Expansion 
Recovery Strain is Also Apparent After Each Spike. 


It is envisioned that future material characterization studies will focus on these new test types 
which extract more fundamental material information, while providing tests which more 
accurately duplicate service environments. 
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3.0 SINGLE CRYSTAL MODELS 


3.1 INTRODUCTION 

The need for accurate constitutive modeling stems from the desire to have accurate life 
prediction models. This goal clearly depends upon understanding the failure process under the 
complex loading and temperature regimes found in gas turbines. However the life models all 
depend on the state of stress and/or strain at critical component locations as input parameters. 
Thus a major part of the life analysis rests on the results of a finite element analysis of the 
component with an accurate material model to give accurate estimates of the stress and strain 
at critical sites in the component. 

The objective of this research is to develop a useful, accurate constitutive model to predict the 
stress and inelastic strain response of the single crystal alloy Rene N6 under the load and 
temperature conditions found in high performance gas turbines, and implement this model in a 
finite element code that can be used for structural life prediction. 

This section contains a review of the relevant material characteristics of Rene N6, the modeling 
approach (discussed in detail in Appendix B), and a few examples of results for Rene N6. The 
section ends with discussion of the additional work that is needed to satisfactorily complete the 
above objective. 

3.2 REVIEW OF RELEVANT MATERIAL CHARACTERISTICS 

The elastic and inelastic properties, and deformation mechanisms characteristic of cubic 
crystals including Rene N6 are much different than traditional polycrystalline alloys as described 
in Appendix A. Both elastic and inelastic response exhibits cubic symmetry. The mechanical 
response depends on the orientation of the applied loads and the properties in tension and 
compression are different when the temperature is below and above about one half of the 
absolute melting temperature. As temperature increases the anisotropy decreases and the 
material is close to isotropic as the temperature approaches the melting temperature. 

The elastic properties change significantly with orientation and Poisson's ratio can even change 
sign. The inelastic properties also exhibit strong orientation dependence with inelastic effects 
like creep, stress relaxation and strain recovery, and the response depends on the rate of the 
applied loads. Even at the same temperature and strain rate the tensile response can vary 
widely with orientation. The shapes of the tensile curves are all different and the ultimate stress 
can vary by 100% or more. 

The active deformation mechanisms by which inelastic deformation occurs fall into three 
general classifications: slip, diffusion, and rotation (see Appendices A3 and A4). Slip results 
from shear stresses propagating dislocations on the crystallographic planes, and active slip 
systems result in the tension-compression and orientation properties. Diffusion processes 
increase with increasing temperature resulting in climb of dislocations and the formation of new 
gamma-gamma prime structures or rafts. Rotation of the crystallographic planes can occur 
during slip to accommodate the deformation and boundary conditions simultaneously when the 
crystal is loaded in a single slip orientation. This type of effect can occur at the blade 
attachment site. Rotation changes the local shear stresses and active slip systems and can 
lead to very large deformations. 


NAS A/CR— 2006-2 13378 


14 



3.3 MODELING APPROACH 


A unified constitutive equation similar to that discussed in more detail in Appendix C is used to 
model Rene N6. Since at least some slip is present at all temperatures the basis of the 
approach is to model the slip in each slip plane and each slip direction. Cubic crystals have 
four octahedral planes with two sets of three slip directions, and three cube planes with one set 
of slip directions on each plane resulting in thirty possible slip systems (see Appendix A4). The 
number of active slip systems depends on the temperature and orientation of the applied load. 
The slip system equations in Appendix A4 are written in the principal coordinates of the cubic 
crystal, so all external loads and boundary conditions must be transformed to the principal 
coordinates of the crystal lattice before using the slip equations. 

There are three slip system constitutive models, one for each of the three sets of slip directions. 
The tension-compression asymmetry results from extended dislocations in the first octahedral 
slip system and the shear stresses resulting in the asymmetry are part of the slip system model. 
The first octahedral slip system model and cube slip system model together give the orientation 
effects. Creep at high temperature and low stress produce slip in the second octahedral slip 
system and this effect is modeled with third the constitutive equation. 

Three sets of state variable evolution equations are used with three slip system models (see 
Appendix A5). There is a back stress variable to model balance between strain hardening 
dynamic recovery observed in tensile and cyclic loading. Temperature effects are included 
through a scalar drag stress that is developed to give the initial hardness (resistance to inelastic 
flow) and long thermal softening resulting from diffusion. The maximum value of the back 
stress (or extent of dislocation pileups) strongly depends on temperature and the local stress on 
the slip plane in the slip direction. The third state variable is used to characterize the maximum 
back stress in creep and high strain rate loading. 

A systematic method has been developed to determine the constitutive parameters from 
tensile, creep and fatigue data (see Appendix A6). The slip equations (Appendix A4) are used 
to determine the local stresses and slip rates in each active slip system from the observed test 
data. These local values are used with the inelastic flow (slip) model and evolution equations in 
a specific order to evaluate and verify the parameters at each step in the process. The 
parameter evaluation method has been automated through the use of EXCEL spread sheets 
and FORTRAN codes. 

3.4 RECOMMENDATIONS FOR FUTURE WORK 

The results to date have been very encouraging, and models for the tensile and tertiary creep 
response are in place. However these results are for one orientation, <001 >, and one 
temperature, 2000°F. To achieve the objective, it is necessary to achieve the following results: 

1 .) Experimentally evaluate the response in other orientations to determine the anisotropy 
parameters in the model. It is necessary to have a complete working model to evaluate the 
response in the blade attachment site. Thus, it is necessary to have tensile, compressive 
and creep data in the <01 1 > and <1 1 1> orientations and different strain rates to determine 
all the anisotropy parameters. Further it would be helpful to have some tensile and creep 
data in other orientations to verify the model. 
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2. ) Extend the mode! to other isothermal temperatures above and below 2000°F. This is the 

first step in developing a model for simultaneous mechanical and thermal cycling. Above 
2000°F the response will be more dependent on temperature than stress. The slip 
orientation effects will become less important but orientation dependent raft microstructures 
will develop from preferred paths of diffusion. Below 2000°F slip becomes more important 
and the slip models may have some different features above and below 2000°F. 

3. ) Extend the model to include thermal history effects and over temperature excursion effects. 

The first step is to use the results for modeling the response at each of the isothermal 
temperatures described in 2) above. This will include evaluating the thermomechanical 
response by interpolation. However our experience has shown that this approach will not 
be totally adequate since the interaction between low temperature slip and high temperature 
climb can produce extra hardening that is not seen in isothermal testing. These studies will 
require the some special thermomechanical experiments to isolate the extra hardening 
effects. 

4. ) Finite element implementation of the model, which is the last step necessary to establish a 

whole blade structural model. An earlier version of the model has been implemented in a 
twenty node brick. The new constitutive model will be implemented in the twenty node brick 
and used with finite element model of the temperature history distribution to establish the 
structural response of a blade and it attachment. The modeling approach will be 
experimentally verified with a simple structure. 
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4.0 GVIPS MODEL FOR DS MATERIAL 


4.1 MODEL BASICS 

A detailed evaluation of the GVIPS model and the COMPARE code is given in Appendix B. As 
described in the appendix, the model has been under development for several years and has 
been applied successfully by Arnold, et. al.. In addition to the basic model, a code called 
COMPARE has been developed which performs nonlinear optimization to determine the 
constants required by the model. As shown in the appendix the model is quite nonlinear, and it 
is convenient to have an automated tool with the capability to perform this non-linear 
optimization. As opposed to the development for the single crystal model given in Section 3.0, 
GVIPS is based on a potential function theory, and consequently, the anisotropic predictions 
are determined by parametric representations in the transverse and longitudinal directions. 
Since the anisotropic behavior predicted by a single crystal model is different than the behavior 
that would be observed for a DS material, the GVIPS model was thought to be appropriate for 
the DS Rene’ 142 data that was tested during this program. Unfortunately, since the testing 
was entirely in the longitudinal direction, the anisotropic predictions of the model were not 
evaluated during this phase of the contract. 

4.2 MODEL CONSIDERATIONS 

As discussed in Sections 2.0 and 3.0, the turbine blade materials being considered in this 
program exhibit a substantial amount of tertiary creep. Prior to this program, GVIPS, like most 
constitutive models did not have a method to predict tertiary creep. During the course of this 
program, Dr. Arnold of NASA Lewis Research developed the capabiltiy to predict tertiary creep 
in the GVIPS model. 

4.3 MODEL PREDICTIONS 

Figures 4-1 through 4-8 show a correlation obtained from COMPARE of some 
Rene’ 1 42 data tested at 1 600°F. Note that these are the results from four tensile and four 
creep tests that were used in getting the correlation. COMPARE has a limitation that only eight 
experiments can be evaluated at a time. As shown, with the exception of the data for the high 
strain rate (Figure 4-8) the correlation is pretty good. The data was analyzed using the version 
of the theory that does not have tertiary creep. Consequently, the data and the predictions 
deviate at longer times (an attempt was made to cull out the tertiary creep data before the 
COMPARE analyses were made; as shown in the figures, some tertiary creep was still included 
in the analyses). 
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Figure 4-1. Comparison Of Predicted Versus Data For Creep Test 57, 
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Figure 4-7. Comparison Of Predicted Versus Data For Tensile Test 5, 
Rene’ 142 At 1600°F - IE-6 In/In/Sec. 
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APPENDIX A - SINGLE CRYSTAL MATERIALS 


A.1 INTRODUCTION 


Nickel base single crystal alloys are used in rocket and gas turbine engines because their high 
temperature properties are superior to those of polycrystalline nickel base superalloys. In high 
temperature applications, grain boundaries in polycrystalline alloys provide passages for 
diffusion and oxidation. Eliminating grain boundaries, and the grain boundary strengthening 
elements, produces materials with superior high temperature fatigue and creep properties 
compared to conventional superalloys. However the absence of grains make the alloys 
anisotropic or orientation dependent. 

This appendix contains a discussion of the elastic and inelastic anisotropic properties of single 
crystal alloys and the details of the development of the constitutive equations that are used to 
predict their response. The appendix begins with an overview of the elastic properties and 
models and the tensor structure of the equations. Sections A.3 and A.4 contain a review of the 
material properties and the three general classes of deformation mechanisms (slip, diffusion, 
and rotation) that control the mechanical properties. The kinematic equations to determine the 
stresses on the slip planes and the accumulated inelastic strain from the slip in each of the 
active slip planes are given in Section A.5. Section A.6 describes the constitutive equations 
used to model slip and diffusion at the slip plane level, and the method used to determine the 
material parameters is given in Section A.7. 

A.2 MATERIALS WITH ELASTIC CUBIC SYMMETRY 


Single metallic crystals with face centered cubic or body centered cubic lattice structure exhibit 
anisotropic elastic response (Dame and Stouffer 1985, 1988). Materials with cubic symmetry 
can be shown to have only three independent elastic constants. The three independent elastic 
constants are designated the elastic modulus, shear modulus and Poisson's ratio. The elastic 
compliance matrix is 


Y v 

r,z 


yz 


1 — V — V 

~E ~E ~E 

-v 1 -v 

~E ~E ~E 

-V -v 1 

~E ~E ~E 

0 0 0 

0 0 0 

0 0 0 


0 0 
0 0 
0 0 
I 0 

G 

0 I 

G 

0 0 


0 



(A2.1) 


This matrix appears identical to the matrix of elastic constants for isotropic materials. However, 
the elastic constants E, v, and G are all independent, and the components of the matrix are a 
function of orientation. Equation A2.1 is only valid when the loading is in the principal directions 
of the material or parallel to the edges of the fee or bcc cubic lattice. Therefore it is necessary 
to establish a method to transform Equation A2.1 to other orientations. The results will be used 
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to determine the elastic constants from three tensile tests on test samples with arbitrary 
orientation relative to the principal axes of the material, and to determine the compliance matrix 
in other orientations. 

Let us establish the elastic stress strain response of a test specimen with the arbitrary 
orientation e' relative to the crystal axes. The vector e* is oriented in space by two successive 
rotations. First is a positive rotation 6 about the e 3 axis so that the e ; basis in the principal 
direction of the material is transformed into the e' basis. Second is a negative rotation (j) about 
e' 2 transforming e' into e'. The two transformation matrices are 


e l 


COS0 

sin0 

o' 

~ e l 

e 2 

= 

-sin0 

COS0 

0 

e 2 

_ e 3_ 


0 

0 

1 

_ e 3_ 


e l 


COS(j) 

0 

sin 4> 

e l 

e 2 

= 

0 

1 

0 

e 2 

_ e §- 


— sin (}) 

0 

COS<j) 

_ e 3_ 


(A2.2a) 


(BA2.2b) 


The constitutive Equation (A2.1) relative to the principal axes of the crystal can be written in 
matrix notation as 

[y'Hs'K] where [S'] = [R][QlR] — 1 [S][Q] —1 (A2.3) 

The transformation matrix [Q] can be established by expanding the transformation rule defined 
by Equation (A2.3) in a six dimensional vector with the components arranged as 

G \i 

G 22 


CT 13 

°23. 

Evaluating the transformation matrix for a rotation about the x 3 axis as described in Equation 
(A2.2a) gives 



c 2 0 

s 2 0 

0 

2c0s0 

0 

0 

s 2 0 

c 2 0 

0 

-2c0s0 

0 

0 

0 

0 

1 

0 

0 

0 

-c0s0 

c0s0 

0 

c 2 0 -s 2 0 

0 

0 

0 

0 

0 

0 

C0 

-s0 

0 

0 

0 

0 

S0 

c0 
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where notation c and s are used for cosine and sine, respectively. The inverse transformation 
is obtained by letting [q - * = [Q(— 0)] » and it can also be shown that [R][Q][R] * = [Q] ^ so the 

transformed compliance matrix can be expressed as 

[S'] = [Q] -T [S][Q] -1 ( A2 - 4 ) 

The compliance matrix can be determined in the e' basis by substituting [Q] _1 into Equation 
(A2.4). The non-zero components in the compliance matrix are 


Si 1 = ± = I( c ° s4 0 + sm 4 0) + ( ^ - f ) C ° S 2 6 sin2 0 


1 2v'\ 9 9 

cos z 0sin l 

G E ) 


S ]2 =— — = ( — Icos^ Osin^ 0 - — (cos 4 0 + sin 4 0 ) 

iZ E' Ie G J EV / 


S14 = — = (-- 2 - - —l cos 0 sin 9 (cos~ 0 -sin^ 0] 
14 E' Ig E j V > 


1 ( 1 + v^\ 9 9 1 

S44 = — = 8 cos z 0 sin z 0 + — 

44 G' \ E J G 


^ 0sin^ 0 + — (cos^ 0 - sin" 0 


Equation (A2.5) shows when the material is loaded with a tensile stress in the e' 3 = e 3 direction, 

a' = [o 0 a 0 0 0 o] T , the transverse strains are independent of orientation and the test 

specimen remains round. If a tensile specimen is loaded in the direction, 

a' = [a 0 0 0 0 0 0] T , then the diametrical response does depend upon orientation. The 

elastic modulus, E 1 , and shear modulus, G’, vary with orientation and Poisson s ratio and v , 
even changes sign. There is coupling between the shear and normal stress and strain 
components defined by the Lekhnitski coefficients of mutual influence. In any arbitrary 
orientation, e '' the matrix of elastic constants can become fully populated and include shear 

stress-shear strain coupling. The shear stress-shear strain coupling constants are called the 
Chentsov coefficients. Thus with only three independent elastic constants, the elastic response 
of a cubic symmetric material is significantly more complex than the elastic response of an 
isotropic material. 

The axial compliance in the e' direction can be determined so that e'n =Sfic 0 for a tensile 
stress vector a" = [c 0 0 0 0 0 0]^. Following the same steps for the second rotation, 
Equation (A2.2b), the compliance of a tensile specimen in the e' direction can be found from 
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Since there are only two independent material parameters in S' u and S]^, from Equations (A2.5) 
and (A2.6) respectively, 


E 


and 


e; 


tensile compliance measurements from the principal orientation and two other different 
orientations can not be used to determine the three elastic parameters E, G and v. Recalling 
that a round specimen remains round when loaded in one of the three principal directions, a 
tensile test in a principle direction can also be used to measure E and v. A second tensile test 
in an arbitrary orientation then gives the shear modulus G. The accuracy of the parameters 
depends on the reliability of the orientation measurements. 

These results were established for evaluating the elastic constants for materials with cubic 
symmetry. Transforming the stiffness or compliance matrices in structural calculations is 
awkward. It is frequently easier to transform stress or strain from the global coordinates to the 
principal coordinates of the material, use the constitutive equation in the principal coordinates of 
the material, and then transform the result back to the global coordinates. Care should be 
exercised when using both engineering and tensor strain. 


A. 3 INELASTIC SLIP MECHANISMS 


The orientation and tension-compression asymmetry properties can be explained from the 
active slip systems and dislocation mechanisms. There may be up to three families of active 
slip systems depending on the orientation of the applied load and temperature. Milligan and 
Antolovich (1989) reported for PWA 1480 for strain rates other than creep (typically greater 

than to -4 sec -1 ) that: 


• Below about 600°C deformation occurs due to slip on the octahedral planes, [111] planes, in 
the 12 directions similar to the [101] direction. Octahedral slip is shown in Figure A-1 and 
corresponds to Slip Numbers 1 to 12 in Tables A-1 and A-2. 

• Between about 600°C and 850°C deformation occurs simultaneously by octahedral slip and 
cube slip. Cube slip occurs on the three planes similar to [001] in the two directions similar 
to the [011] direction (6 slip systems) as defined in Figure A-2 and Slip Numbers 25 to 30 in 
Tables A-1 and A-2. 

• Above about 850°C the strong effect of climb and cross slip reduces the dependence on 
orientation and the particular slip systems that are active. 
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Tensile Loading 


T 



Figure A-1. Effect Of Tensile And Compressive Stresses On The Core Width 
Of An Extended Dislocation. 



Figure A-2. Definitions Of The Stresses, r. 
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Table A-1 . 


Distribution Of The Resolved Shear Stress In Each Slip System Of An FCC 
Crystal For A Loading Of 100 MPa In The [100], [110], [111], And [321] Directions. 


Slip No . 

[ 100 ] 

[ 110 ] 

[ 111 ] 

[ 321 ] 

1 

408.25 

408.25 

.00 

349.93 

2 

.00 

- 408.25 

.00 

- 174.96 

3 

408.25 

.00 

.00 

174.96 

4 

- 408.25 

.00 

.00 

- 116.64 

5 

- 408.25 

.00 

- 272.17 

- 291.61 

6 

.00 

.00 

- 272.17 

- 174.96 

7 

408.25 

.00 

- 272.17 

.00 

8 

.00 

.00 

.00 

.00 

9 

408.25 

.00 

- 272.17 

.00 

10 

.00 

- 408.25 

- 272.17 

- 349.93 

11 

- 408.25 

- 408.25 

- 272.17 

- 466.57 

12 

- 408.25 

.00 

.00 

- 116.64 

13 

- 235.70 

235.70 

.00 

.00 

14 

471.40 

235.70 

.00 

303.05 

15 

- 235.70 

- 471.40 

.00 

- 303.05 

16 

- 235.70 

.00 

- 314.27 

- 269.37 

17 

- 235.70 

.00 

157.13 

33.67 

18 

471.40 

.00 

157.13 

235.70 

19 

- 235.70 

.00 

157.13 

.00 

20 

471.40 

.00 

- 314.27 

.00 

21 

- 235.70 

.00 

157.13 

.00 

22 

471.40 

235.70 

157.13 

336.72 

23 

- 235.70 

235.70 

157.13 

134.69 

24 

- 235.70 

- 471.40 

- 314.27 

- 471.40 

25 

.00 

353.55 

471.40 

454.57 

26 

.00 

353.55 

.00 

151.52 

27 

.00 

353.55 

471.40 

404.06 

28 

.00 

353.55 

.00 

202.03 

29 

.00 

.00 

471.40 

252.54 

30 

.00 

.00 

.00 

- 50.51 
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Table A-2. Correspondence Between The Slip Number In Table A-1 And Slip Directions On 
The Crystal Planes. 


Slip No. 

Slip Plane 

Slip Direction 

1 

Octahedral Slip a/2<l 10>{ 111} 
[111] 

[10-1] 

2 

[111] 

[0-1 1] 

3 

[111] 

[1-10] 

4 

[-1 1-1] 

[1 0-1] 

5 

[-1 1-1] 

[110] 

6 

[-1 1-1] 

[0 1 1] 

7 

[1-1-1] 

[110] 

8 

[1-1-1] 

[0-1 1] 

9 

[1-1-1] 

[10 1] 

10 

[-1-1 1] 

[Oil] 

11 

[-1-1 1] 

[10 1] 

12 

[-1-1 1] 

[1-10] 


Octahedral Slip a/2<l 12>{ 111} 


13 

[111] 

[-1 2-1] 

14 

[111] 

[2-1 -1] 

15 

[111] 

[-1 -1 2] 

16 

[-1 1-1] 

[12 1] 

17 

[-1 1-1] 

[ 1 -1 -2] 

18 

[-1 1-1] 

[-2-1 1] 

19 

[1-1-1] 

[-1 1-2] 

20 

[1-1-1] 

[2 1 1] 

21 

[1-1-1] 

[-1 -2 1] 

22 

[-1-1 1] 

[-2 1 -1] 

23 

[-1-1 1] 

[ 1 -2-1] 

24 

[-1-1 1] 

[112] 


Cube Slip a/2<110>{ 100} 


25 

[10 0] 

[0 1 1] 

26 

[10 0] 

[0 1-1] 

27 

[0 1 0] 

[10 1] 

28 

[0 1 0] 

[1 0-1] 

29 

[0 0 1] 

[110] 

30 

[0 0 1] 

[-1 1 0] 


These results are consistent with the observations of Miner (1 986, 1 986b) for Rene N4 at 760, 
870, and 982°C. Miner found only cube slip for loading near the [111] orientation, and only 
octahedral slip for loading near the [100] and [110] orientations. 

The number of active slip systems similar to [1 1 1 ][1 01 ] depends on the orientation and 
magnitude of the applied load. The shear stress, t, parallel to any slip system can be calculated 
from the applied tensile stress, a, using Schmid's Law 

i = o cos(a) cos(P) (A3. 1 ) 
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where (3 is the angle between the tensile axis and the slip plane normal vector and a is the 
angle between the tensile axis and the slip direction. It can be shown, for the slip systems 
similar to [1 1 1 ][1 01], that loading in the [001] direction will activate eight slip systems 
simultaneously if the stress is equivalent to the critical resolved shear stress (CRSS). Loading 
in the [110] direction will activate four simultaneous slip systems, whereas only one slip system 
will be activated for loading in the [123] direction. 

Next, slip on the octahedral plane occurs by the disassociation of a total dislocation into two 
partial dislocations. The atoms in the core of the extended dislocation have a different stacking 
sequence. The stacking fault energy controls the core width, and the small core widths have a 
greater potential for cross slip to a cube plane. In an ideal setting slip occurs by adding atoms 
to the core and removing atoms from the core. In practice, extended dislocation movements 
are disrupted by other dislocations, antiphase boundaries, cross slip, and climb at high 
temperatures. 

In single crystals, Lall, Chin and Pope (1979) showed that extended dislocations create the 
tension-compression asymmetry because they behave differently in tension and compression. 
Consider a tensile load that creates a shear stress on an octahedral plane that is not parallel to 
the slip direction. This shear stress can be resolved into two components with one component 
parallel to the slip direction, x , and one component normal to the slip direction, Xj , as shown in 
Figure A-1 . During tensile loading in the [001] direction, for example, the normal component of 
the shear stress, x l5 will increase the potential that atoms are added to the core, increase the 
resistance to cross slip, and increase the yield stress. When the applied stress is reversed, the 
shear stress -x^ will favor adding atoms to the core and reduce the potential that atoms are 
removed from the core. Thus the core width is extended during tensile loading and reduced 
during compressive loading. Next recall that at temperatures below the critical temperature 
there are no stress in the cube slip systems. Thus dislocations that cross slip from the 
octahedral planes to the cube planes during tensile loading in the [001] direction are pinned and 
their movement is restricted. At high temperatures the dislocations are not pinned on the cube 
planes and cross slip increases the potential for slip in the cube slip systems. Thus, the core 
width controls the potential for cross slip, and the response in tension and compression will be 
different. 

The movement of extended dislocations can be influenced by other stresses in the crystal 
lattice. Lall, Chin and Pope (1979) also proposed that the stress on the cube plane, x 2 , parallel 
to x on the octahedral plane will promote cross slip from the octahedral plane to the cube 
plane, see Figure A-2. The potential for the cross slip is increased with increasing 
temperatures. The stress x 2 is particularly important when the loading is near the [111] 
orientation where the stresses on the cube plane may be larger than the octahedral stress. 
Paidar, Pope and Vitek (1984) later proposed that the stress normal to the slip direction on the 
intersecting octahedral plane (X3) also influences the core width. The effect of X3 is similar to 
that of x t on the core width, and increasing values of x 3 increase the potential for cross slip of 
the dislocation. 

Milligan and Antolovich (1989) investigated the tensile and cyclic behavior of PWA 1480 in the 
temperature range from 20°C to 1093°C at strain rates of 0.50 and 50.0% per minute. They 
found two types of dislocation networks at high and low temperatures. At temperatures 
between approximately 20 and 760°C and for loading in the [001] orientation with plastic strains 
up to 0.3 %, the dislocation network was dominated by pairs of extended dislocations shearing 
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the y ' phase. In this temperature range dislocations traveled through the y / as closely spaced 
pairs. In contrast to the shearing mechanism at low temperature, a y ' bypass mechanism was 
observed at temperatures of 982°C and higher. For specimens tested to failure at 982°C there 
were high dislocation densities in the y phase, and that deformation was controlled by looping 
or bypass of the y' particles. Shearing of the y' phase was also observed after a large 
increase in the matrix dislocation density and significant strain hardening. In the intermediate 
temperature range from 760°C to 927°C a transition from shearing to bypass was observed. 

To summarize, deformation of single crystal occurs by slip in the octahedral and cube slip 
systems. The tension-compression asymmetry results from the propagation of partial 
dislocations that depend on four shear stress components, x and %2 in the slip direction, and 
and x 3 that are normal to the slip direction. To use this information in constitutive modeling it is 
necessary to develop the equations for local stresses from the applied stress in a global or 
reference coordinate system. It will also be necessary to determine a representation for the 
inelastic strain resulting in the global coordinate system from the local slip in each slip system. 

A-4 LATTICE ROTATION AND CREEP 

The deformation discussed so far has been limited to small plastic deformations typical of low 
cycle fatigue applications. However, if the deformations are larger than 3 or 4%, such as found 
in creep and tensile loading, rotation of the crystallographic structure may occur during slip. 

The rotation occurs in single slip orientations to simultaneously accommodate the deformation 
and restraints at the boundary. For example, consider a tensile test of the single crystal tensile 
specimen oriented in a single slip orientation. The grips hold the ends of the specimen in a 
fixed alignment. During application of the load, slip will initiate on a number of parallel slip 
planes located close to the maximum shear stress, and slip bands will be observed on the 
surface of the specimen. As deformation progresses, continued slip would not be possible 
without rotation of the slip planes to accommodate the change in length of the specimen. 
Bending will occur near the grips to maintain continuity and alignment of the specimen. During 
elongation the slip planes rotate toward the tensile axis, and the angle between the slip plane 
and tensile axis decreases. As rotation continues a second slip plane will become oriented in a 
position to permit slip. When the stress on this plane reaches the critical resolved shear stress, 
slip will occur in a second slip system. Eventually the slip on the primary and secondary planes 
will become equal and duplex slip will occur (identical slip in both systems). 

The rotation of a specimen in a single slip orientation can be plotted in a stereographic 
projection. Specimen D in Figure A-3a has one active slip system in the [ 101] direction for 
loading in the [123] as direction shown in Table A-1 . The tensile force produces a rotation that 
moves the [101] slip system vector toward the tensile axis at Point D, and the angle X 
decreases in Figure A-3a. However, instead of rotating the stereographic projection toward 
Point D, it is more convenient to leave the stereographic projection fixed and rotate Point D 
toward the [Toi] slip system. Using this approach, when Point D reaches Point 2, two slip 
systems ([1 1 1][l01]] and [1 1 1][01 1]) will have the same stress and duplex slip will result. From 
Point 2 the specimen will rotate along the [001]-[1 1 1] boundary to [1 12], the point midway 
between [Toi] and [Oil] in the [ill] plane (the vector sum of the two slip systems). The 
specimen will remain at this orientation until failure. Figure A-4a shows the directions of 
rotation along the [001 ]- [ 1 11] boundary for a specimen starting in any orientation within the 
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stereographic triangle. Crystals with an initial orientation on the perimeter of the stereographic 
triangle will exhibit duplex slip upon loading, and rotate initially to the poles shown in 
Figure A-4a. 



(a) ( b ) 


Figure A-3. Stereographic Projection Showing The Rotation Of A Tensile Specimen Loaded 
In A Single Slip Orientation. 



[U2] 

* 211 ] 

(o) ( b ) 

Figure A-4. Directions Of Rotation For Crystals Loaded In Single Slip Orientations In The 

Standard [001] Stereographic Triangle For The (A) [1 1 1] [1 1 1] And (B) [1 1 1] [1 12] 
Slip Systems. 

Creep deformation at the higher stress levels, but below the yield stress, can occur by 
dislocation creep and dislocation glide. At lower values of stress, creep bulk diffusion occurs 
(Nabarro-Herring Creep). Single crystal creep is more complicated because the active slip 
system for dislocation creep or dislocation glide can change during creep or be different for 
different alloys. Leverant and Kear (1969, 1970, 1973) and MacKay and Maier (1982) reported 
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for Mar M200 and Mar M247, respectively, that primary creep and secondary creep occur in 
different slip systems. In particular, for temperatures near 750°C: 

• Primary creep occurs by deformation in the slip systems similar to [1 1 1][1 12]. This is the 
second octahedral slip shown in Table A-2 and corresponds to Slip Numbers 13 to 24 in 
Table A-2. 

• Secondary creep occurs by deformation in the slip systems similar [1 1 1 ][1 01 ]. 

However, Hopgood and Martin (1986) reported for SRR99 at 750 and 800°C that: 

• Both primary and secondary creep occurred by deformation in the slip systems similar to 

[111H112]. 

• Deformation in the cube system has not been reported for creep near the [111] orientation. 

Thus, it appears that primary creep occurs in the slip systems similar to [1 1 1][1 12], but the 
active slip system for secondary creep may be an alloy specific property. Lattice rotations for 
the [1 1 1][1 12] slip system are shown in Figure A-4b. The boundary X-Y is defined by the 
Schmid Factor cos(A)cos(0) = 0.427 in Equation A3-1 . For the [1 1 1][1 12] slip system, duplex 

slip only occurs along the [001 ]-Y and X-[Tll] boundaries of the stereographic triangle. 

Nickel base single crystals also exhibit a delay in the onset of creep for small values of stress. 
Pollock and Argon (1991) conducted a study of the creep resistance of CMSX-3 in the [001] 
orientation at temperatures between 800°C and 900°C and loads between 450 MPa and 600 
MPa. There was no instantaneous straining on loading (elastic strains are subtracted out), and 
delays of 1 .6x1 0 4 and 600 seconds were recorded for the tests at 800°C and 850°C, 
respectively. No delay was observed at 900°C, and lower stresses were found to produce 
longer incubation periods in the onset of primary creep. To document the process during the 
incubation period, a creep test was terminated during the delay period at 4.5x1 0^ seconds. A 
typical micrograph shows that dislocations spread out to dislocation free regions from initial 
grown-in dislocation tangles. Primary creep is assumed to be initiated when an active 
deformation system is achieved by the meeting of dislocation loops. 

A. 5. KINEMATIC EQUATIONS FOR CRYSTALLOGRAPHIC SLIP 

It is necessary to develop the equations for local shear stresses and local slip in each slip 
system for use in constitutive modeling (Dame and Stouffer, 1985, 1988). The stress vector t 
acting on some plane with an outward normal vector n is given by 


V 


’°11 

CT 12 

a 13 

’nf 

t2 


°21 

CT 22 

a 23 

n 2 

_ c 3 . 


_ a 31 

CT 32 

a 33_ 

_ n 3 _ 


(A5.1) 


where ay are the stress components in a principal coordinate of the material (parallel to the 

right hand coordinate system in the [100], [010] and [001] orientations) at the point in question. 
The component of the stress, x , in some direction p is then calculated from 
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1 ^ = [ t l *2 t 3 ] ' [PI P2 P3] 


(A5.2) 


where p and n are unit vectors. The local shear stresses can be determined from Equations 
A5.1 and A5.2 for each of the 30 slip numbers given in Table A-1 and shown graphically in 

Figures A-1 and A-2. The result for the primary octahedral slip systems similar to [lll][ 1 01] is 
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The shear stresses in the slip systems similar to the secondary octahedral [1 1 1 ][1 1 2] are given 
by 
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and the stresses in the cube slip systems similar to [001 ][01 1 ] are 
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The final step is to develop a correlation between stress in the slip direction, x, and the stresses 
tj, x 2 , and x 3 that control the core width of a partial dislocation in the octahedral slip system. 

Consider, for example, the stresses in the [1 11] [10 1 ] slip system. Recall that x t is normal to x 
in the slip plane, thus can be identified as x 13 . The stress X 2 is parallel to xj or [ 10 1 ] on an 
intersection cube plane, thus X2 is % 28 - The stress x 3 is normal to x in the intersecting 
octahedral plane; thus x 3 is xjg. The same approach is used to determine the remaining 
correlations given in Table A-3. 


Table A-3. Correlation Among the 12 Shear Stresses r in the Octahedral Slip System 
and the Stress Components x,, r 2 , andx 3 That Control the Core Width of a 
Partial Dislocation. 


T 

T i 

T 2 



1 

13 

28 

16 


2 

14 

26 

20 


3 

15 

30 

24 


4 

16 

28 

13 


5 

17 

29 

19 


6 

18 

25 

22 


7 

19 

29 

17 


8 

20 

26 

14 


9 

21 

27 

23 


10 

22 

25 

18 


11 

23 

27 

21 


12 

24 

30 

15 



The inelastic strain in the principal coordinates of the material will be developed next from the 
slip in each of the slip systems. The local slip is equivalent to the local engineering shear 
strain. Let y represent the local engineering shear strain in the right hand coordinate system 
defined by the unit normal vector n, the unit slip direction vector, p, and s = n x p where s is the 
vector parallel to the slip systems similar to [1 1 1][1 12]. Applying the second order tensor 
transformation rule gives the global tensor strain as a function of the local slip in the form 
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Calculation of the inelastic strain, e- 1 , in the principal coordinates of the material due to the slip 
in the 12 primary octahedral slip directions gives 
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and the strain in the secondary octahedral [1 1 1 ][1 12] are given by 
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The strain, e??, for the cube slip system can be written as 
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As a check of the calculations, notice that the inelastic strains given by Equations A5.7 to A5.9 
are symmetric. 

A.6 MODELING OF SINGLE CRYSTAL ALLOYS 

The next step in the process is to develop a constitutive model to relate the slip on the slip 
plane to the shear stress on the slip plane in the slip direction following the methods developed 
by Sheh and Stouffer, 1988 and 1990. Let us begin by reviewing the ideas that will be used to 
develop a model: 

Only octahedral slip is observed at temperatures below the critical temperature, about 700°C, 
for high strain rate loading. At temperatures above the critical temperature, deformation occurs 
by both octahedral and cube slip. Cube slip dominates near the [ill] orientation and 
octahedral slip dominates near the [001] and [01 1] orientations. 

Slip, or plastic flow, occurs when the shear stress on the slip plane in the slip direction, t, 
exceeds the critical shear stress. The tension-compression anisotropy occurs because the 
stresses normal to the slip direction on the slip plane and intersecting slip plane, t| and %2 
respectively, extend and contract the core width of partial dislocations. Small core widths 
increase the potential for cross slip. The shear stress on the cube plane parallel to an active 
slip system, x 3 , also promotes cross slip at higher temperatures. 

The small strain tensile response show the presence of strain recovery, thus a back stress 
variable should be included as part of the model. The variation in recovery with orientation 
shows that the effect of back stress will be different in the octahedral and cube slip systems. 

The flat tensile curves show that there must be a balance between strain hardening and 
dynamic recovery. The strain rate sensitivity occurs from a change in the balance between the 
hardening rate and recovery rate. 

The cyclic hardening suggests the need for a drag stress variable to model cyclic hardening 
that is different from strain hardening. However, note that some two-phase superalloys are 
nearly stable in cyclic testing, so only the initial value of the drag stress may be required. 

Crystal lattice rotation occurs above a few percent strain in single slip system orientations. 
Crystal lattice rotation effects are expected to be present in tensile and creep response beyond 
3 to 5% strain. 
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The inelastic or plastic strain or strain rate is obtained by summing the slip or slip rate of the 
individual slip systems. The strain in the principal coordinate of the material for the octahedral 
and cube slip systems is given in Section A.5. 

The constitutive model is based on a unified strain approach where the total strain is 
decomposed into elastic, inelastic and thermal components. The coordinate system is 
assumed to be the principal coordinates of the material. For fee crystals the principal directions 
are along the cubic edges or in the [100], [010] and [001] directions. Thus the load or strain in a 
structure must first be transformed to the principal coordinates of the material before using the 
equations developed in this section. The elastic strain in the principal directions is given by 
Equation A2.1 . Recall for cubic symmetry that the elastic modulus, shear modulus and 
Poisson's ratio are independent. A method for determining the elastic constants was presented 
in Section A.2. The thermal strain in the principal coordinate system is isotropic; therefore the 
thermal strain components can be determined from the coefficient of thermal expansion. 
Equations A5.7, A5.8, and A5.9 can be used to determine the total inelastic strain rate as 


C< ~/3 

£ ij - £ ij + £ ij + e ij 


(A6.1) 


for the octahedral and cube slip systems acting simultaneously. This approach will require 
three separate flow equations, one for each type of slip system. Recall that partial dislocations 
are present in octahedral slip but not in cube slip. The parameters must be determined such 
that the correct flow equation will be operative under the correct conditions as summarized 
above. The accumulated inelastic strain is determined by integrating the inelastic strain rate. 
The modeling task therefore reduces to determining the flow equations for the octahedral and 
cube slip systems, incorporating the tension-compression asymmetry, and evaluating the 
material parameters. 

Several of the polycrystalline materials have rather similar chemical compositions to the y -y' 
alloys and exhibit many similar time dependent features. Thus it is reasonable to use the same 
type of equation for the slip rate in the octahedral and cube slip systems. However, recall that 
the stresses in all three-slip systems will be required to model the core width effect. Defining 
the following indicial notation for the three slip systems: 

k = 1 defines slip systems similar to [1 1 1 ][1 01 ] and a = l, 4 and p = 1, 3 
k = 2 defines slip systems similar to [1 1 1][1 12] and a = 1,4 and p = 1,3 
k = 3 defines slip systems similar to [001 ][01 1 ] and a = 1,3 and p = i,2 


Defining the inelastic flow equation in the above notation gives 



= D k exp 
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n kl 





(A6.2) 


where and are the drag stress and back stress, respectively, that are associated with 
the slip rate in each of the slip systems. The shear stresses in the active slip systems are 
denoted by and the difference between the shear stress and back stress, ' s 
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net stress required to produce slip. The parameters D k and A k are scale factors that will be 
used to separate the high and low strain rate response, and n kl controls the strain rate 
sensitivity. 

The back stress evolution equations discussed in appendix C included strain hardening and 
dynamic recovery terms. Sheh (1988) found a back stress model that satisfied the saturation 
condition and was more flexible than those in appendix C for modeling the response of the 
three types of slip systems. Writing his result in the above notation gives 
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where F k> G^, n k3 and £2 k at are material parameters. The inelastic recovery is primarily 
modeled by the strain rate term. If the Gk term is small relative to the Fk term, the changes in 

the value of the back stress will be small during elastic loading and unloading. If there is a 
large change in the back stress during the inelastic strain recovery, then contribution of the 
second term is relatively small in the [001] orientation. The exponents in Equation A6.3 give 
more flexibility to model the amount of recovery. The parameter x 0 = 1 MPa or KSI is included 

to normalize units of stress. The saturated back stress, £2 s k at , is assumed to be a constant for 

tensile and cyclic loading. 

The drag stress is used to model the initial hardness through Z 0 and the growth in hardness 

through an exponential evolution equation. Drag stress was chosen to modify the tension 
compression anisotropy because the yield stress is a measure of the initial resistance to slip. 
Using the slip system notation the cyclic growth in drag stress can be rewritten as 

H k P =H k2 -(H kl -H k2 )exp^-m k w“ P j (A6.4) 


where = } r^d y , is the accumulated inelastic work in each slip system. The 
parameters Hki and Hk2 are the initial and saturated values of the drag stress, respectively. 
Equation A6.4 is a measure of the resistance to slip in each slip system due to the initial 
hardness and the development of a dislocation network. Since the tension-compression 
asymmetry occurs in the slip systems similar to [1 1 1 ][1 01 ], only the needs to be modified. 

The simplest approach is to assume a linear relation between z®^and x®^, x®^, and x®^, i.e., 1st 


z °cp 


= h“ ( 3 +Vix“ ( 3 +v 2 


P ap 


+v 3 tf 


(A6.5) 


The material parameters V- 1 , V2 and V3 must be determined to give the required asymmetry. 
From their definitions, x?^ and x?^ are defined to extend or constrict the core width and the 

1 ^ ctB 

ability of the dislocation to cross slip from the octahedral to cube plane. Thus when x, p and 
xf are positive the resistance to slip, Z®^, is increased, but the resistance to slip is decreased 
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when xf and x“ ^ are negative. Thus the resistance to inelastic flow depends on both the 
magnitude and sign of x^ and x^. Recalling that x^ is the thermally aided contribution of the 
cross slip, and the sign of x“ p is not relevant since the potential to cross slip is the same in 
tension and compression. The parameters V-| , V2 and V3 become less important with 
increasing temperature above the critical temperature. At very high temperatures Zj ^ — > Hj ^ 
and the orientation dependence disappears. The stresses x“^, x“^ and x“^ are defined in 
Table A-3. 

A.7 MATERIAL PARAMETERS AND CALCULATED RESULTS 

The following general assumptions and limitations are made to guide the development of the 
material parameters: 

• The material constants are determined under isothermal conditions. Extension to 
thermomechanical response has not been attempted for single crystal alloys. 

• The material parameter evaluation procedure is for strain rate controlled tensile and cyclic 
data. 

• Single crystal response depends on strain rate and temperature. At high strain rate below 
about 600°C deformation occurs in the 12 slip systems similar to [1 1 1 ][1 01 ], between about 
600°C and 850°C deformation occurs simultaneously by octahedral slip and cube slip (6 slip 
systems similar to [001][011]), and above 850°C orientation dependence becomes 
negligible. 

• The octahedral slip constants are determined for tests were there is no cube slip. For 
example, data from tests in the [001] and [01 1] orientations can be used at any temperature 
to determine the octahedral constants. Conversely, the cube slip constants were 
determined from tests in the [111] orientation. The local stresses are given in Table A-1 for 
a 1000 pound tensile load in the [001], [011], and [111] orientations. 

Table A-1 shows that tensile loading in the [001] orientation will produce equal slip rates in eight 
slip systems. The shear stress in the slip plane producing slip is 0.408 of the applied tensile 
load. Likewise, there are assumed to be eight equal slip rates, y . that add up to the inelastic 

strain rate, e 1 , in the [001] orientation. The slip system indices are dropped since there are 
eight equal shear stresses and slip rates in the [001] orientation. In this case Equation A5.7 
can be used to show that y = (l 2 /l 6 j 6 )e l = 0.306£' after noting that stresses corresponding to 

the Y21 » y\i' Y 42 and Y43 s,i P systems are negative (see Table A-1). Similar results can be 

established for the [011] (or [111]) orientation since the are four (or three) equal shear stresses 
and slip rates. Thus, the local shear stresses and local slip rates for tensile, cyclic and creep 
tests in the [001], [011] and [111] orientations can be obtained from experimental data. 

Since the equations are the same in both octahedral and cube slip systems, the discussion of 
determining parameters will be limited to the octahedral system as an example. This system 
includes the tension-compression asymmetry terms. The method of determining parameters is 
divided into the following three sequential steps for the evaluation of the: 
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• flow equation parameters 

• back stress parameters 

• drag stress parameters. 

The order of the steps must be preserved since the information flow is sequential. 

In the flow Equation A6.2 the parameter D 0 is a somewhat arbitrary scale factor that correlates 
with the maximum inelastic strain rate. Typically it is a factor 10 4 sec.’ 1 greater than the 
maximum expected strain rate in a particular application. For example, D 0 = 1 sec.’ 1 for creep, 
D 0 = 10 4 sec.’ 1 for tensile and fatigue response for strain rates up to 1 .0 sec.’ 1 , and 
D 0 = 10 6 sec. 1 or greater for wave propagation and high rate loading applications. 

Since n controls the strain rate sensitivity, let us begin with the parameters in the flow equation 
and experimental tensile response data. At tensile data saturation, (the strain at which the 
stress-strain curve becomes flat) the inelastic strain rate is equal to the total strain rate, 
which is constant in a constant strain rate tensile test. At saturation the back stress is at the 
maximum value, and the drag stress is near the initial value Z 0 and remains almost 

constant. Thus for saturated, isothermal, tensile loading Equation A6.2 becomes 


Y L - A ex p| 


-A, 


■'lo 




(A7.1) 


without the slip system numbers since they are all equal. The last term, x -Q/|r -n| , in Equation 
A7.1 is equal to plus one. Using this result, Equation A6.2 can be inverted and written in the 
form 


X l / n u = 

1 

In 

fr u) 


. A > 


^ A J_ 


li 


Ao 

T; — ^1/n 


(A7.2) 


which can be used to find n, £2 m and Z 0 assuming that d { is assigned a value as described 
above. Equation A6.5 can be rewritten as 

t s xf" = H„ + V,Ti s + v 2 |t 2s | + VjT* +£2 S X S (A7.3) 


after using Equation A6.4 in Equation A7.2 with k = 1 . Recall the initial value of the drag stress 
in the octahedral slip system is Hi ■) , and Q. s is one of the eight equal values of the saturated 

back stress in the [001] orientation. Assuming a value forn^ (the strain rate sensitivity 
exponent), Equation A7.3 becomes a linear equation with five unknown parameters, Hi i , Vi , 

V2, V3, and £» s . At least five sets of data pairs ( X^ 2 " , t s ) are required to evaluate these 

parameters in a least squares routine. These data can come from tension and compression 
tests at different strain rates in the [001] and [011] orientations. An iteration procedure can be 
used to determine the parameters. To summarize, start by assuming a value for ni 1 
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( 0 . 1 > nu > 3 . 0 ) and calculate H-ji, V-|, V 3 , V 3 , and Q sat from the data. Iterate until the 
saturated back stress falls in an acceptable range as described next. 

The maximum value of the back stress can vary widely in different types of materials. The 
value during tensile loading must be positive (opposite to the tensile stress) and less than the 
saturated tensile stress, that is 0 < £2 m < o s . In an extensively warm worked material that has 

developed subgrains and raised the yield strength considerably, the back stress is quite small 
when compared to the drag stress (Miller, 1975). In undeformed large grain materials the 
maximum value of the back stress can easily be greater than one half of the saturated tensile 
stress. Milligan and Antolovich (1989) determined values for the saturated back stress as large 
as 0.85 g s in PWA 1480, a single crystal alloy with large precipitates. Large values of back 

stress are also verified by cyclic loading if reverse plastic flow begins before the applied stress 
reaches zero on unloading. To summarize, the following assumptions can be use if no other 
data is available: 

Q. m =0.0 to 0.3a s for small grain materials ^ ^ 

£2 m =0.5 to 0.7 o s for 1 arge grain materials 

It is very difficult to determine from only tensile data; therefore it may be necessary to use 
an approximation. If cyclic data is available then other methods can be used to determine £2 m . 

The parameters in the back stress equation will be determined next. If only tensile data is 
available it is reasonable to assume the maximum value of the back stress can be 
approximated by Equation A7.4, and the values of Z 0 and n have been determined as described 

above. In this case a back stress history can be determined from the tensile data by using the 
inverted flow equation, Equation A7.2. The back stress history can then be used to determine 
the back stress rate as a function of time, 

Q(t) = 40 i7S7r[ H l 1 + V l*l w + V 2 h(t)| + V 3 T 3 W] (A7.5) 

X v (t) 


The back stress rate can be obtained by numerical differentiation of the back stress history. 

This is reasonable for a computer acquired data file with hundreds of data points, but it not easy 
for hand digitized data. 

The constants g i and n 12 in the back stress Equation A6.3 can be determined from the initial 
value of the back stress rate during the initial “elastic" loading phase of the tensile history. 

Since the initial back stress rate is proportional to the stress rate or strain rate, 


Qi(0) = G, 


f r \ 


V X o ) 


ir,(0) 


(A7.6) 


the initial value of the back stress rate can be correlated to the local stress rate. The 
parameters G x and « 12 are determined from the log linear slope of the line. The parameters 

n -13 and F-j for the inelastic term can be determined from the back stress rate in a tensile test. 
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The back stress rate is then determined numerically using the post yield uniaxial test data. 
Rearranging the Equation A7.5 for k = 1 and x 0 = 1 MPa or KSI gives 

= Fj| yp 3 (A7.7) 


The equation can be used with tensile data and/or recovery to determine ni 3 and F-| . A log- 
linear curve of the data can be used to determine the parameters n-| 3 and F-| . 

The next step is to develop a method to evaluate the material parameters for the drag stress 
equation. First, recall that the parameters in the flow and back stress equations have been 
evaluated. The inverted flow Equation A7.3 can be solved for the drag stress history 

H x {t) = T X 1/2n - Vjr, - V 2 |t 2 | - V3T3 -Q X (A7.8) 

for log term tensile loading or cyclic loading. The parameter X 1//2n is the mapped inelastic 
strain rate defined in Equation A7.2. 

The parameter m t that controls the rate of hardening or softening can be estimated by 
estimating the total accumulated inelastic strain or work required to produce the cyclically stable 
response. As the material response approaches the long-term stable response the exponential 
function in Equation A6.4 must approach zero and H x (t) = H i2 . Assuming a stable response 
occurs when the drag stress is within one percent of saturation, Equation A6.4 becomes 

— ^2- = 0.01 = exp(-/n 1 lT, / j (A7.9) 

where w/ is the accumulated inelastic work required to reach the cyclically stable state. The 
value of m, can then be estimated from the accumulated inelastic work at saturation. The 

accumulated inelastic strain |e‘| can be substituted for W 1 to obtain a slightly different 
representation. 
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APPENDIX B - EVALUATION OF THE GVIPS CONSTITUTIVE MODEL AND COMPARE 
CODE 

B.1 INTRODUCTION 

Generalized Viscoplasticity models with a Potential Structure (GVIPS) have been developed to 
predict the mechanical response of high temperature materials with anisotropic and inelastic 
properties under complex theromechanical loading environments, (Saleeb and Wilt, 1993, and 
Arnold and Saleeb, 1994). Further these models have been implemented in the computer code 
COMPARE to evaluate the material parameters from files of experimental data using nonlinear 
optimization techniques, (Saleeb et. a!.). The purpose of the current study is to evaluate the 
possible application of the GVIPS model for the directionally solidified ahoy R142 at high 
temperature. 


B.2 THE GVIPS MODEL 

The GVIPS model is comprised of coupled system of tensorial rate equations which give the 
evolution of the state variables as a function of the external loading. The flow equation for the 

Inelastic strain rate £ 1 is given by 



and the evolution law for the back stress state variable cc is given by 


(B.2.1) 


a = 


Z.+ 


h' 


h(l + 2j3) x - 


a®a 


(h£'-R„G m n) 


j v 


(B.2. 2) 


where h- H/G p , h' = -fi H/G p+ ' and Z m = M *. The other parameters in the inelastic flow 
and back stress equations are given by 


1! 

cP 
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P" 
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i 

(B.2. 3) 

2 k 0 " 


G = J- T (a\.M:(a) 

(B.2. 4) 



T = M : (cr— a j 

(B.2. 5) 

n = M:a 

(B.2. 6) 

M = P-£Q--gR 

(3.2.7) 


In the application for R142 alloy the isotropic hardening parameter k o is held constant. The 

anisotropy of the material is introduced through Equation B.2. 7. The two parameters £ and C, 
are defined as 
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(B.2.8) 


and 


1 = 


t] z -l 

2 

V 


with 



K, 



with 



(B.2.9) 


In the above equations, the material parameters and k, are the longitudinal and transverse 
shear strengths and, Y t and Y t are the longitudinal and transverse strengths in tension. The 
fourth order tensors P, Q and R are defined as 

P = /-i(<S®s) (B2 - 10) 

Q m =\(d,jS m +D„S lk + D ll S n +D,Sik)-W :J D l , (B-2.11) 

R = 3(d®£>)-(d®5+5® d)+j(s®5) (B.2.12) 

The tensor D is used to define the principal direction of the material as a function of the 
direction cosines, and 8 is Kronecker delta. The GVIPS model has nine material parameters 
( k a , n, n, m, (3, R, H, t, and Q that must be determined for any application. The inelastic strain 
rate must be integrated and added to the elastic strain to find the total strain. 


B.3 EVALUATION OF A UNIAXIAL GVIPS EQUATION 


The first step in determining the uniaxial model is to evaluate the tensor M for loading parallel to 
the principal axis of the material. In this case the material behavior will be pseudo-isotropic and 
only the tensor P in Equation B.2.7 is required. Let A be a second order tensor that represents 
the three normal components of stress or strain. The relationship between a three dimensional 
physical space and a six dimensional mathematical space is the given as 


a tl 0 
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a 33 


= - a i 


a 33 0 0 0; 


(B.3.1) 


Loading along the principal axis of the material will give the six-dimensional vector as 
D = [ 1 O' 0 0 0 o] and the Kronecker delta vector is written as <5 = [l 1 1 0 0 0], 

so the Equation B.2.10 becomes 
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after eliminating the higher order tensor components which are all equal to zero. Equation 
3.3.2 is the representation for an isotropic material. Finally note that Equation B.3.2 is singular, 
so Z m =M -1 in the back stress equation can not be evaluated. This appears to be a 

fundamental inconsistency in the potential function formulation. The authors make the 
assumption that 


Z m =M~ l 


1 G 0 
0 1 0 
0 0 1 


(6.3.3) 


for use In three dimensional applications. 

The next step is to use Equations B.3.2 in Equations B.2.3 to B.2.7 to obtain a model for 
uniaxial loading. In the case we assume the uniaxial stress produces a non-zero inelastic strain 
rate in the three coordinate directions. This implies that there are three non-zero back stress 
coordinates, but the symmetry of the loading implies that e 22 = £ 33 = -e n /2 and 

a, 2 = a 33 =-a u /2. Define s’ , a and a as the uniaxial values of the inelastic strain rate, 
stress and back stress, respectively, then the overstress vector becomes 

|( 7 -aJ= [a -a, -a/2, -a/2] (B.3.4) 


Equations B.2.5 and B.2.6 then become 


r = i(cr-3a/2 
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(B.3.5) 

The scalar functions F and G in Equations B.2.3 and B.2.4, respectively, can be written as 


F = -L(a- 3 a/ 2 f -l 

jK q 

g =A(3 

3k- 

The inelastic flow equation for uniaxial loading is now given as 
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and the back stress equation can be rewritten as 

1 


a =• 


{1 + 2/3) 


f u 3 

e' -aR„G n 


G 


P 


(B.3.6) 

(B.3.7) 


(B.3.8) 


(B.3.9) 


after simplifying the results. 

These equations have been implemented In a simple FORTRAN computer code and integrated 
using a generalized trapezoidal meihod with small time steps. The results match the results 
from COMPARE using constants that were determined by COMPARE. The test cases included 
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tensile loading at constant strain rate, creep under constant stress, and a complicated multiple 
step creep loading. 

B.4 EVALUATION OF THE MATERIAL PARAMETERS 

In order to understand the roll of the material parameters in the GVIPS constitutive model, the 
equations were evaluated under specific conditions. This produced another approach to 
evaluating several of the parameters and helped to generate a single set of parameters that 
could be used for both tensile and creep loading. This result also provides guidance in 
selecting the initial values used in the COMPARE code. 

Saturation . Saturation is defined at the condition when the stress and inelastic strain rate 
become constant in tensile or creep loading. This corresponds to the minimum creep rate 
during creep loading or the ultimate tensile stress during constant strain rate loading. In both 
cases the back stress in Equation B.3.8 is constant and the back stress rate is zero. This 
condition can be used to simplify Equations B.3.8 and B 3.9. 

Consider first the inelastic flow equation. Since the term in the ( ) brackets is much greater 
than zero at saturation, Equation C.3.8 can be approximated as 

3^(3 K;)e‘ s ={a s -3a s /2f n+l (B.4.1) 

where the subscript s corresponds the value of the corresponding variables at saturation. 
Define the following parameters 

D = (3k ~ ) and p= — — (B.4.2) 

so the inelastic flow equation becomes 

(A,+ t=(<T s -3a: 5 /2) (B.4.3) 

This can be written an a logarithmic linear representation in the form 

p\ogD 0 +p\ogs‘ s — log (o' 5 -3a s /2) (BAA) 


which can be solved for D 0 and p using the method of least squares and the saturation data if 
the saturation values of the back stress are known. 


The back stress evolution rate will be zero at saturation, thus Equation B.3.9 can be written as 
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after using Equation B.3.7 and rearranging. Once again define two additional combinations of 
constitutive constants as 
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so the back stress equation at saturation can be rewritten as 



(B.4. 6) 
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Applying the method of least squares to Equation 8.4.7' in the form 

log e' s = log P a +q log (3a' v /2) (B.4.8) 

gives another linear representation that the creep and tensile saturation data must satisfy. 

Back Stress History . Assume a strategy can be found to estimate the values of the back stress 
at saturation, a s , so that D 0 and p can be found from the saturation data. Defining 

X{f)={p o e‘ (t)) P as a transformed value of the inelastic strain rate at any time in a deformation 
history, Equation 8.4.3 can be solved for the back stress history 

X{t)] (B.4.9) 

3 


upon using the stress and inelastic strain rate history from any test data file. The back-stress 
rate history, a(t), can then be estimated by fitting a spline to the back stress history. 

Rearranging Equation 3.3.9 and using Equation 8.4.6 gives 
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which is another representation that the parameters must satisfy using data from a tensile or 
creep test. Define Q 0 as 

Q 0 = (1 + 2/3)/# (B.4.11) 


and observe Equation C.4.10 can be written as 
log Q„ + 8 log 
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(B.4.12) 


Equation B.4.12 can be solved for Q 0 and (3 using the method of least squares after estimating 
a value for Ko. 

A Solution Strategy . At this point there are six parameter combinations (D 0 , p, P 0 , q, Q 0 , P) to 
determine seven constitutive parameters (jx, Ko, n, H, P, R«, m) using ihe method of least 
squares. Further it is necessary to obtain an estimate for the saturated value of the back 
stress, a s before any of the equations can be used. The saturated value of the back stress is 
typically less than about 40% of the stress at saturation. The value of saturated back 
stress, a s , physically corresponds to the balance between strain hardening and dynamic 

recovery, and depends on the grain size or mean distance between obstacles that pin 
dislocations. The back stress is usually negligible in alloys with very small grains, and near the 
upper limit for large grain alioys like R142. A good initial guess for R142 would be 30% to 35% 
of the saturated stress. Thus one approach is to estimate the value of a s , and iterate on the 

estimate until the calculated values of the saturated stress and inelastic strain rates optimized. 
This method has been successfully used with other constitutive models, Stouffer and Dame, 1996. 

A second approach was developed to estimate the saturated back stress as part of the effort to 
verify the above equations to evaluate the material parameters. In this case several R142 
creep and tensile test files were loaded into the COMPARE code to determine the material 


NASA/CR— 2006-2 13378 


51 



parameters. The parameters were then used in the uniaxial code to evaluate the back stress 
history. It was found that the saturated back stress, a s , increases iineariy from about 25% of 

the saturated creep stress to about 35% of the saturated back stress from the tensile tests. 
This saturated back stress model was used with the above equations to determine the material 
constants as described above. The calculated parameters were very close to the parameters 
determined form COMPARE, and the calculated model results were identical. This exercise 
confirms the above method to determine the parameters and provides insight into a method to 
model the saturated back stress for R142 at other temperatures. 

Finally the calculated results were determined by assuming a small value for Ko In a separate 
study it was shown that the calculated response is essentially the same for any small value of 
Kq Thus Ko can almost be taken as a universal constant. 

B.5 A MODEL FOR TERTIARY CREEP 

In a separate study S. Arnold (1 998) developed a model for tertiary creep that has been 
implemented into the uniaxial code and evaluated. The parameters in the model appear to 
have sufficient flexibility to model the tertiary creep of R142. 


B.6 OVERALL EVALUATION OF THE GVIPS/COMPARE SYSTEM 


In general, it was found that several changes had to be made to the COMPARE program to 
effectively use the program. These changes mainly embodied changing the user interface to 
facilitate entering experimental data so that multiple runs could be made without manually re- 
entering the data files. Once these changes had been made the program was much easier to 
run. However, COMPARE places several demands on the format of the experimental data. At 
first, these requirements were not understood since the code is still under development and a 
user manual was not available. These requirements include defining strain and time points at 
each point that a load change is made, and including the load up strains or stresses in the file 
(even though they are not always included in the measured data such as in a creep test). Once 
all the requirements were understood, it was relatively easy to run COMPARE and to generate 
meaningful results (and results that could be interpreted consistently with the uniaxial theory 
developed above). 


Additionally, it should be remarked that the GVIPS theory is based on potential theory and not 
on metallurgical principals. As described in Appendix A, the single crystal theory relies heavily 
on such physical principals. It will be interesting to see how well the GVIPS approach will work 
on a total data set that incorporates several temperatures and several deformation 
mechanisms. 


NAS A/CR— 2006-2 13378 


52 



B. REFERENCES 


Arnold, S.M. and Saleeb, A.F. (1994) "On the Thermodynamic Framework of 

Generalized Coupled Thermoelastic-Viscoplastic Damage Modeling", Int, Journal of Plasticity, 

V 10.3, pp 263-278. 

Saleeb, A.F. and Wilt, T.E. (1993) "Analysis of the Anisotropic Viscoplastic Damage Response of 
Composite Laminates" Int. J Numerical Methods Eng. V.36, pp. 1629-1660. 

Stouffer, D.L. and Dame, T.L. (1996). Inelastic deformation of metals: Models, mechanical properties, and 
metallurgy, John Wiley and Sons Inc., New York. 


NASA/CR— 2006-2 13378 


53 




APPENDIX C - MODELING APPROACH 


C.1 CHARACTERISTICS OF MODERN CONSTITUTIVE EQUATIONS 

The state variable approach resulted from a synthesis of modeling, physical metallurgy and 
experimental results of many researchers and practitioners in the last two decades. The 
development was motivated by three general conclusions about the state of classical modeling. 
First, the yield surface approach to plasticity is not totally adequate. It is excellent for predicting 
the initial yield of metals, but the subsequent yield surfaces that evolve during a general 
multiaxial deformation history are not easy to model or experimentally quantify. The plastic 
strain increment vector is not always normal to the yield surface nor is it parallel to the 
deviatoric stress vector. Second, the addition of a creep strain to account for time effects is not 
fully adequate. The coupling between plastic and creep strain can be significant even in the 
traditional "constant stress" creep test where it is known that the initial loading rate effects the 
subsequent creep response. Time is not a good variable for constitutive modeling because it is 
not objective. Third, the results of physical metallurgy are largely ignored. An enormous body 
of fundamental knowledge about metals has been developed during the last half century. 
Information about the physical basis of deformation and the associated variables should be 
useful for modeling. Knowledge of the active deformation mechanism under different loading 
conditions will provide insight about the macroscopic response and clues for the choice of 
variables. A prior knowledge of the deformation mechanisms will guide experimental programs 
and the selection of materials for specific applications. 

This section contains a discussion of the variables that are used in the state variable modeling 
and the structure of the equations that form the philosophy of the approach. The results from 
physical metallurgy and the classic approaches are examined to identify some specific choices 
for the state variables and the effect of these choices on the mathematical structure of the 
equations. 

Characterizing the properties of dislocations and their reactions with grain boundaries, 
precipitates, and other dislocations could take many variables. Alternatively, implicit internal 
state variables are defined as macroscopic averages of events associated with the 
development of the microstructure. They characterize average properties associated with 
dislocation propagation, and their interaction with precipitates, grain boundaries, and other 
dislocations. In addition to the state variables there are external variables that are 
macroscopically defined or observed, and controlled, measured or calculated in an 
experimental program. These variables must be correlated with the interna! state variables of 
the material during deformation. 

Plastic and creep strain variables are not used. Instead a single, unified, inelastic strain 
variable is used to characterize all the inelastic effects. The kinematic equation used to 

determine the total strain tensor, gT is the sum of the elastic, gE, inelastic, gl, and thermal 

strain tensors, gTH. Mathematically, the kinematic equation is represented by 


ej =gy +e!j+g™ (C.1.1) 

Equation C.1 .1 avoids many of the issues associated with defining and characterizing plasticity 
separately and developing a rule for creep-fatigue interaction. The inelastic strain can be 
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uniquely determined experimentally from the total measured strain since the elastic and thermal 
strains are almost entirely unaffected by inelastic deformation. Time is avoided as a variable. 
Time dependent variables will be characterized by their rates such as the inelastic strain rate or 
state variable evolution rate. The inelastic strain rate, for example, can be integrated to obtain 
the accumulated inelastic strain at any point in the deformation history without ever using time 
explicitly as a state variable. 

One objective of state variable constitutive modeling is to develop a correlation between the 
externa! variables and the internal state variables. The existence of this correlation permits the 
inelastic strain rate of a metal at any point in the deformation history to be determined from the 
current values of the external and internal variables. The inelastic strain rate at a later point in 
the history is determined by later values of the mechanical and thermal loads, and the later 
values of the state variables. Thus, the state variables must be continuously updated to 
account for the changes in the microstructure. The relationship between the current 
microstructure and previous deformation history is included in the model by updating the state 
variables; that is, the "history dependence", is included through the evolution of the state 
variables! The entire previous history of the deformation is characterized by the current state of 
the microstructure. 

Mathematically, the state variable approach can be characterized by three types of equations: 
a kinematic equation, a kinetic equation, and a set of state variable evolution equations. The 
kinematic equation is used with the integral of the inelastic strain rate, that is, 


£ X 

& ij 


= e y + J^ 

o 


ydt + eif 


(C.1.2) 


where t is some time in the deformation history. The total strain and inelastic strain can be 
large or small, but the elastic strain is assumed to be small so a linear elastic stress-strain law 
can be used. The kinetic equation, or inelastic fiow law, is formally written as 


£ ij — ^ij kl ’ X> ’ ’ 


•) 


(C.1.3) 


where s' is the inelastic strain rate tensor that must be integrated with respect to time to obtain 
the accumulated inelastic strain. The stress, a,, , temperature, T, and other external and 
interna! state variables . are expected to be present. The state variables are used as 
macroscopic measures of the characteristics of the microstructure. The evolution equations to 
update the current values of the state variables are formally given as 

tk=4k(oij.E|j.T.?i.-) <C - 1 - 4) 

and must also be integrated with respect to time to determine their current values. At this point 
in the development the exact definition or number of state variables is undefined and needs to 
be established from the principles of physical metallurgy. Equations (C.1 .2) to (C.1 .4) establish 
the essential mathematical structure for all state variable models. 

The next step is to identify some potential state variables. The state variables must be selected 
to describe the critical microstructural effects in anticipation that the predicted macroscopic 
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inelastic response will be accurate to a broad class of loading conditions and temperatures. 
Thus iet us review the critical metallurgical and experimental observations that can be used to 
define the state variables and external variables. These important concepts are: 

a. ) Inelastic deformation occurs at low temperature due to slip. The slip occurs on the slip 

planes closest to the planes with the maximum shear stress. 

b. ) Inelastic deformation occurs at high temperature due to planar slip and dislocation 

climb, a diffusion controlled mechanism. 

c. ) The accumulated inelastic strain does not produce a significant change in the inelastic 

volume. 

d. ) The accumulated inelastic strain and strain rate appear to be independent of hydrostatic 

stress except for some alloys under extreme hydrostatic loading conditions. 

e. ) Back stress (a resistance to slip) results from the interaction of dislocations with other 

dislocations, grain boundaries, precipitates, and other barriers. The net stress 
producing slip or inelastic strain is the over stress or the difference between the shear 
stress and back stress. The back stress is in the direction opposite to the local shear 
stress in uniaxial loading and is therefore expected to be orientation dependent in three 
dimensional loading. 

f. ) The initial hardness of a metal is determined from the resistance to slip that is created 

by obstacles such as grain boundaries, precipitates, and interstitial atoms in the initial 
microstructure. These effects are isotropic and are modeled by the initial value of the 
drag stress. 

g. ) Hardening is an increase in the resistance to deformation. It results primarily from the 

interaction of dislocations with other dislocations, precipitates and grain boundaries. 

h. ) Dynamic recovery is the annihilation of dislocations during deformation that reduces the 

effective rate of hardening. It was shown that temperature and strain rate dependence 
results from the balance between the hardening and recovery rates. 

i. ) Recovery is achieved by the formation of cells and subgrains. Cells and subgrains are 

orientation independent and monotonically approach a saturated state under uniform 
loading conditions. 

j. ) Creep results primarily from the temperature dependent mechanisms o? diffusion, 

dislocation glide, climb and rafting in single crystal superalloys. Creep in tension and 
compression are different due to the effect of lattice expansion or compression on the 
diffusion rate. 

k. ) Static recovery results from the interaction stress between the dislocations themselves. 

The higher temperatures increase dislocation mobility and promote diffusion. Static 
recovery lowers the effective steady state hardening during creep. The rates of dynamic 
and static thermal recovery are much different. 
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The net effect of hardening and softening can be separated mathematically into two groups by 
their orientation dependent or orientation independent characteristics. This establishes the 
need for a minimum of two types of hardening state variables. A back stress tensor is 
introduced to model the hardening and recovery effects associated with the back stress effect. 
In addition, a scalar drag stress is used to model the development of hardening or softening 
associated with isotropic effects such as point defects, precipitates, cells and subgrains and 
grains. Other scalar staie variables are introduced to mode! the net effect of hardening and 
static thermal recovery during creep. These variables will be introduced as necessary in the 
development of the different models. 

The back stress effect has been incorporated in state variable models by many workers. In all 
cases, the kinetic equation for the inelastic strain rate, Equation C.1 .3, has the form 

where the deviatoric stress and deviatoric back stress are used to preserve the condition of no 
inelastic volume changes and the function X depends on the other macroscopic and state 
variables. Equation (C.1. 5) includes the effects of back stress by permitting reverse inelastic 
flow after unloading when the over stress, defined as the difference between the applied stress 
and the back stress, becomes negative. The tensorial property of the back stress also allows 
for orientation dependence that is consistent with the dislocation pile-ups discussed above. If 
the applied stress direction is changed, the direction of shear stress in the crystallographic slip 
system will also change, and slip would occur in another siip plane in a different direction 
producing different dislocation interactions. Thus the components of the back stress tensor will 
also change. Finally, the initial value of the back stress is usually assumed to be zero because 
dislocation interactions do not develop until after slip has started. 

The drag stress variable, Z, is introduced to model the initial hardness, and the subsequent 
hardening or softening that develops during cycling. The initial value of drag stress, Z(0)= Z 0 , 

is used to characterize the initial resistance to slip that is present in the undeformed 
microstructure. The evolution of the drag stress is associated with the interaction of 
dislocations with point (orientation independent) defects and the growth of ceils and subgrains. 
Since the interaction of dislocations with precipitates, interstitial atoms, cells and subgrains is 
largely orientation independent, the drag stress is assumed to be a scalar. The drag stress is 
assumed to monotonically increase (for hardening) or decrease (for softening) to a saturateo 
state in constant amplitude cycling. 
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